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Abstract. There has been much success in describing the hmiting spatial fluctuations 
of growth models in the Kardar-Parisi-Zhang (KPZ) universality class. A proper 
rescaling of time should introduce a non-trivial temporal dimension to these limiting 
fluctuations. In one-dimension, the KPZ class has the dynamical scaling exponent 
z = 3/2, that means one should find a universal space-time limiting process under the 
scaling of time as iT, space like t^l^X and fluctuations like t^l'^ as i — ^ oo. 

In this paper we provide evidence for this belief. We prove that under certain 
hypotheses, growth models display temporal slow decorrelation. That is to say that in 
the scalings above, the limiting spatial process for times t T and t T+t^ are identical, for 
any v < 1. The hypotheses are known to be satisfied for certain last passage percolation 
models, the polynuclear growth model, and the totally / partially asymmetric simple 
exclusion process. Using slow decorrelation we may extend known fluctuation limit 
results to space-time regions where correlation functions are unknown. 

The approach we develop requires the minimal expected hypotheses for slow decor- 
relation to hold and provides a simple and intuitive proof which applied to a wide 
variety of models. 



1. Introduction 

Kardar, Parisi and Zhang (KPZ) [26] proposed on physical grounds that a wide variety 
of irreversible stochastically growing interfaces should be governed by a single stochastic 
PDE (with two model dependent parameters D,X ^ 0). Namely, let x h- )■ h{x,t) G M 
be the height function at time t and position x G M'', then the KPZ equation is 



dh{x,t) r^\ui \ \\\-!Ui m2 



dt 



DAh{x,t) + X\Vh{x,t)\^ + r]{x,t), 



where ri{x, t) is a local noise term modeled by space-time white noise. Since then, it has 
been of significant interest to make mathematical sense of this SPDE (which is ill-posed 
due to the non-linearity) and to find the solutions for large growth time t. 

Significant progress has been made towards understanding this equation in the one- 
dimensional d = 1 case. Specifically, it is believed that the dynamical scaling exponent 
is 2 = 3/2. This should mean that for any growth model (also polymer models) in 
the same universality class as the KPZ equation (i.e., the KPZ universality class), after 
centering h by its asymptotic value h{v) := \imt-^oo \h{vt,t) and rescaling^ 

h{vTt + X{Ttfl\ Tt) - Tt h{v + X(Tt)-V^) 
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the limit of ht{X,T) should exist (as t — )■ oo) and be independent of v. Moreover, the 
limit should, regardless of microscopic differences in the original models, converge to 
the same space-time processjn. 

Most of the rigorous worlcl done in studying the statistics associated with this fixed 
point have dealt with the spatial process (obtained as the asymptotic statistics of 
ht{X,T = 1) as a process in X, as t — oo) and not on how the spatial process evolves 
with T. The exact form of these statistics depend only on the type of initial geometry 
of the growth process (e.g., the Airy^ process for non-random flat geometries and Airy2 
process for wedge geometries; see the review [20]). 

Computations of exact statistics require a level of solvability and thus have only been 
proved in the context of certain solvable discrete growth models or polymer models in 
the KPZ universality class. The partially/totally asymmetric simple exclusion process 
(P/TASEP), last passage percolation (LPP) with exponential or geometric weights, the 
corner growth model, and polynuclear growth (PNG) model constitute those models 
for which rigorous spatial fluctuation results have been proved. Recently, progress was 
made on analyzing the solution of the KPZ equation itself [HIZllISllSl] , though this still 
relied on the approximation of the KPZ equation by a solvable discrete model. 

The slow decorrelation phenomenon provides one of the strongest pieces of evidence 
that the above scaling is correct. Indeed, slow decorrelation means that ht{X,T) — 
ht{X, T + 1^~^) converges to zero in probability for any u < 1. Fix m times of the form 
Tt + aif^ (for G M and < i < m). Then, as long as u < 1, the height function 
fluctuations, scaled by t^^^ and considered in a spatial scale of t^/^, will be asympotically 
(as t — )■ oo) the same as those at time Tt. 

Specifically, we introduce a generalized LPP model which encompasses several KPZ 
class models. Then we give sufficient conditions under which such LPP models dis- 
play slow decorrelation. These conditions (the existence of a limit shape and one-point 
fluctuation result) are very elementary and hold for all the solvable models already 
mentioned, and are believed to hold for all KPZ class models. The proof that slow 
decorrelation follows from these two conditions is very simple - it relies on the superad- 
ditivity property of LPP and on the simple observation that if Xt > Yt and both Xt and 
Yt converge in law to the same random variable, then Xt — Yt converges in probability 
to zero (see Lemma [2. 3p . 

Previously, the slow decorrelation phenomenon was proved for the PNG model [TU] . 
Therein the proof is based on very sharp estimates known in the literature only for the 
PNG. Apart from the PNG, the only other model for which slow decorrelation has been 
proved is TASEP under the assumption of stationary initial distribution [4J. 

Besides being of conceptual interest, the slow decorrelation phenomenon is an impor- 
tant technical tool that allows one to, for instance: (a) easily translate limit process 
results between different related observables (e.g., total current, height function repre- 
sentation, particle positions in TASEP; see [3]), and more importantly, (b) prove limit 
theorems beyond the situations where the correlation functions are known [Sl fT^I^ (see 
Section [XT!) . A further application is in extending known process limit results to prove 
similar results for more general initial conditions / boundary conditions [T5] . 

^As explained in the forthcoming paper |17) . this space-time process is expected to be a non-trivial 
renormalization fixed point for the whole KPZ universality class. See also [1S1[571[1S| for previous 
discussion of space-time scalings in the physics literature. 

■^Results for joint distributions at different times before taking t to infinity have been derived, first 
in the problem of tagged particle in the TASEP (35], and then in more general models [5[|T^. However, 
the different times are restricted to lie in an interval of width 0(t^^^). 
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Outline. In Section [2] we introduce the general framework for LPP models in which we 
prove a set of criteria for slow decorrelation (Theorem 12. ip . In the rest of the paper, we 
apply Theorem 12.11 to various models in the KPZ class, which can be related in some 
way with a LPP model: the corner growth model, point to point and point to line LPP 
models, TASEP, PASEP (which requires a slightly different argument since it cannot 
be directly mapped to a LPP problem) and PNG models. Finally we note extensions 
of the theorem to first passage percolation and directed polymers, provided that (as 
conjectured) the same criteria are satisfied. 

Acknowledgments. The authors wish to thank Jinho Baik for early discussions about 
this and related problems. I. Corwin wishes to thank the organizers of the "Random 
Maps and Graphs on Surfaces" conference at the Institut Henri Poincare, as much of 
this work was done during that stay. Travel to that conference was provided through 
the PIRE grant OISE-07-30136 for which thanks goes to Charles Newman and Gerard 
Ben Arous for arranging for this funding. I. Corwin is funded by the NSF Graduate 
Research Fellowship. S. Peche would like to thank Herve Guiol for useful discussions on 
TASEP and her work is partially supported by the Agence Nationale de la Recherche 
grant ANR-08-BLAN-0311-01. The authors are very grateful to the anonymous referee 
for careful reading and a number of constructive remarks. 

2. A SUFFICIENT CONDITION FOR SLOW DECORRELATION 

In this section we consider a general class of last passage percolation models (or 
equivalently growth models). Given the existence of a law of large numbers (LLN) 
and central limit theorem (CLT) for last passage time (or for the associated height 
function), we prove that such models display slow decorrelation along their specific, 
model dependent, "characteristic" directions. 

We consider growth models in M'^^^ for > 1 which may be lattice based or driven 
by Poisson point processes. We define a directed LPP model to be an almost surely 
sigma- finite random non- negative measure /i on W^^^. For example we could take fi to 
be a collection of delta masses at every point of Z'^+^ with weights given by random 
variables (which need not be independent or identically distributed). Alternatively we 
could have a Poisson point process such as in the LPP realization of the PNG model. 
We will focus on a statistic we call the directed half-line to point last passage time. We 
choose to study this since, by specifying different distributions on the random measure 
H one can recover statistics for a variety of KPZ class models. In order to define this 
passage time we introduce the half-line 

nC = {p : Pi = P2 = ■ ■ ■ = Pd+i < 0}, 

where pi is the i coordinate of the point p. 

It is convenient for us to define a second coordinate system which we call the space- 
time coordinate system as follows: Let R be the rotation matrix which takes HC to 
{p Pi < 0,p2 = ■ ■ ■ = Pd+i = 0}. Then the space-time coordinate system is 
applied to the standard basis. The line {p : Pi = P2 = ■ ■ ■ = Pd+i} (which contains HC) 
is the inverse image of {p : p2 = ■ ■ ■ = Pd+i = 0} and we call it the t-axis (for "time"), 
see Figure [1] for an illustration. The other space-time axes are labeled xi through Xd 
(these are considered to be "space" axes). Call a curve tt in M'^^^ a directed path if 
7 = Rir is a function of t and is 1-Lipschitz. Two points are called "time-like" if they 
can be connected by such a path. Otherwise they are called "space-like" . 
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Figure 1. The black line is half-line T-LC and the space and time axis 
are label. A directed path vr from HC to the point p is shown. 



To a directed path we assign a passage time 



T(7r) = fi{7r) 



which is the measure, under the random measure /i of the curve it. Now we define the 



where we understand the supremum as being over all directed paths starting from 
the half-line and going to p. One may also consider point to point last passage time 
between p and q which we write as Lp-p{p,q). This is the special case of /i = on 
{x : X - p eR'^+^\R'i+'^}. 

In Section [3] we show how, by specifying the random measure fi differently, this model 
encompasses a wide variety of LPP models and related processes (such as TASEP and 
PNG). Just to illustrate though, take d = 1 and let fi be composed of only delta masses 
at points p in with mass Wp exponentially distributed with rate 1. Then Lncip) is 
the last passage time for the usual LPP in a corner (or equivalently the corner growth 
model considered in Section [3?T|) . We present our result in this more general framework 
to allow for non-lattice models such as the PNG model. 

We can now state a result showing that slow decorrelation occurs in any model which 
can be phrased in terms of this type of last passage percolation model provided both a 
LLN and a CLT hold. 

Theorem 2.1. Fix a last passage model in dimension d+1 with d > 1 by specifying the 
distributions of the random variables which make up the environment. Consider a point 
p G W^~^^ and a time-like direction u G M!^"^^. // there exist constants (depending on 
p, u, and the model): i^c (^i^d i-p-p non-negative; ^uc-,lvv ^ (0, 1); v G iS^^luchw)'! 
distributions D, D' ; and scaling constants c%Ci cvv such that 



last passage time from the half-line TiC to a point p as 



Ludp) 



sup T(7r) 

tt:HC-^P 




as t goes to infinity, 




Cn^t^nc 
+ f'u) - t^'t 



as t goes to infinity, 




as t goes to infinity. 
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then we have slow decorrelation of the half-line to point last passage time at tp, in the 
direction u and with scaling exponent v , which is to say that for all M > 0, 

lim ¥{\L'Hc{tv + f'u) - Lnc{tp) - f'irrl > Mt^^'') = 0. (2) 

t—^oo 

Remark 2.2. There are many generalizations of this result whose proofs are very 
similar. For instance the fixed (macroscopic) point p and the fixed direction u can, in 
fact, vary with t as long as they converge as t — t- oo. One may also think of the random 
LPP measure n (and the associated probability space fl) as depending on t. Thus for 
each t the LPP environment is given by fit defined on the space flf The probability 
P will therefore also depend on t, however an inspection of the proof below shows that 
the whole theorem still holds with P replaced by F^. 

Proof of Theorem \2.1\ Recall the super-additivity property: 

L-Hcitp + > L-Hc{tp) + Lvv{tp, tp + f'u), 

which holds provided the last passage times are defined on the same probability space. 
This follows from the fact that, by restricting the set of paths which contribute to 
Lnc{tp + f^u) to only those which go through the point tp, one can only decrease the 
last passage time. The following lemma plays a central role in our proof. 

Lemma 2.3 (Lemma 4.1 of [6]). Consider two sequences of random variables {Xn} 
and {Xn} such that for each n, Xn and Xn are defined on the same probability space 
VLn- If Xn > Xn and Xn =^ D as well as Xn ^ D then Xn — Xn converges to zero in 
probability. Conversely if Xn =^ D and Xn — Xn converges to zero in probability then 
Xn =^ D as well. 

From now on, we assume that the different last passage times Ly^ci-) and Lp-p(-) are 
realized on the same probability space. Also, by absorbing the constants c-^£ and cj}-p 
into the distributions, we may fix them to be equal to one. Using super-additivity we 
may write 

Lndtp + fu) = Lnc{tp) + Lvv{tp, tp + eu) + X^, (3) 
where > is a (compensator) random variable. Rewriting the above equation in 
terms of the random variables Xi{t)-, X2{t) and Xsi^) and dividing by t'^'^^ we are left 
with 

X2{t) = xi{t) + x^^ity^^^--^^" + Xit-^«^. 

By assumption on z/, z/7p-p — < and hence we know that the distribution of X2{t) 
and separately of Xi(^) + X3(^)^'''^^^ '''^^ converge to the same distribution D. However, 
since Xft"''""^ is always non-negative, we also know that X2{t) > Xi(t) + X3{t)t^^''^''^~'^'^^ . 
Therefore, by Lemma 12.31 their difference, Xtt~'^'^^, converges to zero in probability. 
Thus X2it) — Xiit) converges to zero in probability. Since 

Luc (tp + ru)- Lnc (tp) - rivv 
X2{t) - Xi{t) = 

the theorem immediately follows. □ 
3. Slow decorrelation in KPZ growth models: 

EXAMPLES AND APPLICATIONS 

The aim of this section is to make a non-exhaustive review of the possible fields of 
applications of Theorem 12. 1[ We introduce a few standard models and explain, briefly, 
how they fit into the framework of half-line to point LPP and what the consequences 
of Theorem 12.11 are for these models. 
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Figure 2. Corner growth height function. The are the random, 
exponentially distributed times it takes to fill an outer corner. The black 
line is the height function at time t = Wi^i. There are two outer corners 
at that time. 

3.1. Corner growth model. We choose to first develop in detail the implications of 
slow decorrelation for a simple KPZ growth process. This process is known as the corner 
growth model and is related to both LPP and TASEP . 

Consider a set At in with initial condition Aq = and evolving under the 
following dynamics: from each outer corner of At a [0, 1) x [0, l)-box is fill at rate 
one (i.e., after exponentially distributed waiting time of mean 1). See Figure [2] for an 
illustration of this growth rule where the model has been rotated by 7r/4. One can 
record the evolution of the growing interface dAt in terms of the random variable 

L{x, y) ■= inf {t > 0| (x - i, y - I) ^ A} for (x, y) G Z^. 

This random variable is well known to define a last passage time in a related directed 
percolation model, as we now recall. Let Wij be the waiting time for the outer corner 
(i, j) to be filled, once it appears. A path vr from (1, 1) to {x,y) is called directed if it 
moves either up or to the right along lattice edges from (1, 1) to {x,y). To each such 
path 71, one associates a passage time T(7r) = j)eiT'^id- Then 

L{x,y)= niax T{tt), 

i.e., L is the last passage times from (1, 1) to {x,y). 

Alternatively one can keep track of At in terms of a height function h defined by the 
relationship 

{h{x-y,t)>x + y} = {L{x,y)<t} (a;,y)GZ^, (4) 
(together with linear interpolation for non- integer values of x — y). Note that for given 
t, h{X,t) = \X\ for |X| large enough. Thus the corner growth process is equivalent to 
the stochastic evolution of a height function h{X,t) with h{X,0) = \X\ and growing 
according to the rule that local valleys (\/) are replaced by local hills (/\) at rate 
one. We will speak mostly about the height function, but when it comes to computing 
and proving theorems it is often easier to deal with the last passage picture. 

3.1.1. LLN and CLT. Analogous to the LLN for sums of i.i.d. random variables, there 
exists an almost sure limit shape for the height function of this growth model: 

Hv) := hm ^^""^'^^ - + ^ 



*-s-oo t forf^(— 1,1) 
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If we more generally consider the height function arising from a LPP model with 
an ergodic measure (e.g., iid random lattice weights), then super-additivity and the 
Kingman's ergodic theorem implies the existence of a (possibly infinity) limit growth 
evolution 

Since LPP is a variation problem, the limiting profile is also given by the solution to 
a variational problem which, when translated into the limiting height function means 
that h satisfies a Hamilton- Jacobi PDE dxh = f{dxh) for a model dependent flux 
function /. Such PDEs may have multiple solutions, and h corresponds to the unique 
weak solutions subject to entropy conditions. Such PDEs can be solved via the method 
of characteristics [18]. Characteristics are lines of slope /' along which initial data 
for h is transported. In our present case if we set p = |(1 — dxh) then p satisfies the 
Burgers equation drp = dx{p{^ — p)) and the characteristic lines are of constant velocity 
emanating out of the origin. 

It is the fluctuations around this macroscopic profile which are believed to be univer- 
sal. Johansson ^23j proved that asymptotic one-point fluctuations are given by 

where -Fgue is the GUE Tracy- Widom distribution defined in [3^. Unlike the traditional 
CLT the fluctuations here are in the order of t^^^ and the limiting distribution is not 
Gaussian. 

Likewise we may consider the fluctuations at multiple spatial locations by fixing 

X(r) =t;t + r(2(l-i;2))V3t2/3^ 

1 -^ 7)2 n _ 7,2^2/3 

Hir, s) = i±^t + r .(2(1 - v')y/H^^^ + (r^ - s) ^' e'\ 

Here H{t, 0) = t h{X{T)/t) + o(l) and H{t, s) — H{t, 0) measures the fluctuations with 
respect to the limit shape behavior. Then, in the large time limit, the joint-distributions 
of the fluctuations are governed by the so-called Airyg process, denoted by A2- This 
process was introduced by Prahofer and Spohn [32j in the context of the PNG model 
(see also [21]). A complete definition of the Airy2 process is recalled in [15]. More 
precisely, it holds that 

\nnF I []{h{X{T,),t) > H{T,,Sk)} ) = P ( f]{A2{rk) < s,} ) , (6) 



^fe=i / \fc=i 



where m > 1, ti < T2 < ■ ■ ■ < and si, . . . ,Sm are real numbers. Of course, ([5]) is the 
special case of ([S]) for m = 1 and ri = 0. 

3.1.2. Slow decorrelation in the comer growth model. We now consider how fluctuations 
in the height function are carried through time. For instance, if the fluctuation of the 
height function above the origin is known at time t (large) for how long can we expect 
to see this fluctuation persist (in the t^^^ scale)? The answer is non-trivial and given 
by applying Theorem 12.11 there exists a single direction in space-time along which the 
height function fluctuations are carried over time scales of order t^, while for all other 
directions only at space-time distances of order t^/^. Indeed, given a fixed velocity 
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V G ( — 1,1), any exponent z/ < 1, and any real number 9, 

{\\h{v{t + et''),t + et") - {t + er)h{v)] - [h{vt,t) -th{v)] \ > Mt^/^) = o, (7) 



t— >oo 



for any M > 0. Thus, the height fluctuations at time t at position vt and at time and 
t + Of^ at position v{t + Of^) differs only of o{t^^^). These fixed velocity space-time lines 
are the characteristic of the Burgers equation above. 

Thus, the right space-time scaling limit to consider is that given in equation ([1]), with 
h and h now taken to be the TASEP height function and asymptotic shape, rather 
than that of the KPZ equation. As noted below equation ([1]), the value of the velocity 
V should not affect the law of the limiting space-time process. As evidence for this, 
equation dH]) shows that we encounter the Airy2 process as a scaling limit regardless of 
the value of v G (—1, 1). This amounts to saying that the fixed T marginals of the full 
space-time limit process of equation ([T]) are independent of v. 

3.1.3. Implications of slow decorrelation. Up to now only the "spatial-like behavior" of 
the space-time process ([T]), i.e., the process in the variable x for fixed /3 (which one can 
set to 1 w.l.o.g.) has been obtained, while the process in the variable /3 remains to be 
unraveled. 

A consequence of ([7]) is that if we look at the fluctuations at two moments of time 
t and t' with \t' — t\ with u < 1, it corresponds to taking /3 = 1 + (9(t''~^) in the 
r.h.s. of ([T]). Then in the t — )■ oo limit, the limit process is identical to the process for 
fixed P = 1. So, if we determine the limit process for any space-time cut such that in 
the t — )■ oo limit, /3 — )■ 1, then, thanks to slow decorrelation, one can extend the result 
to any other space-time cut with the same property. In the following we refer to this 
property as the process limit extends to general 1 + dimensional space-time directions, 
meaning that we have 1 dimension with spatial-like behavior and dimensions in the 
orthogonal direction. 

As indicated in the Introduction, slow decorrelation also allows for instance (a) to 
translate the limit of different related observables and (b) to extend results on fluctuation 
statistics to space-time regions where correlation functions are unknown. We illustrate 
these features in the context of the corner growth model. For simplicity, we consider 
the case where v = 0. Fix m > 1, i/ G [0, 1), real numbers ti < T2 < ■ ■ ■ < Tm and 
si, . . . , Sm- Then set 

-- [l{t + etn-r2-^^''t^/^\, (8) 



/3 



i{T,9,s) = (t + + (s - r2)22/¥ 

(a) We first show that one can recover ([6]) from an analoguous statement in the 
corresponding LPP model using (jl]) and slow decorrelation. We start from a result 
in LPP. Consider the fixed y = t/4 slice of space-time. This is obtained by setting 
Okf = rfe2^/3^2/3 for ^j^jpj^ 

x{t, e) = \t + t2}IH^I\ y{r, 9) = \t, £(r, e,s) = t + t2^'H'"'' + {s - T^)2'/h^'\ (9) 
Using the Schur process [H], it is proven [T5|l22] that for Okf" = Tk2'^^H'^/^, k = 1, . . . ,m, 

Jhn P I f]{L{xiTk,9k), \t) < i{Tk,Sk)} ) = P I f]{A2{Tk) <Sk}]. (10) 
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Figure 3. Assume that the black dots are 0{T'^) for some u < 1 away 
from the line y = t/4. Then, the fluctuations of the passage time at the 
locations of the black dots are, on the T^/^ scale, the same as those of 
their projection along the critical direction to the line y = t/A, the white 
dots. 



To get the result on the height function using (jlj), one would need to make the 
choice: Okt" = -(s^ - t^)2'^/H^/^. Then we have 

x(r, 6) - y{T, 6) = T2'^h^^\ x(r, 6) + y{T, 6) = - {s - t^)2-^'H^/\ £(r, 9, s) = t. 



Thus, to obtain ([6]) (for v = 0) from (|TOi) it is actually sufficient to project (x, y) in ( 1TT|1 
on the line y = t/A along the characteristic line passing through {x,y), see Figure |3] for 
an illustration. One finds that this projection gives the scaling (|9]) but with r replaced 
by some f = r(H-o(l)) — j-rast— i-oo. The reason is that the characteristics for r 7^ 
have slope slightly different from 0. Finally, slow decorrelation (Theorem 12. see also 
d?])) and the union bound imply then that 



(b) The results for and (ITOjl are derived by using the knowledge of (determinantal) 
correlation functions. The techniques used for these models are however restricted to 
space-like paths (in the best case, see |H]), i.e., for sequences of points (x^, yk)k such that 
Xk+i — Xk >0 and yk+i — yk^O (which can not be connected by directed paths). Now, 
choose in ([8]) Okf^ = O^t^ — (s^ — r|)2^/^t^/^ for some real 6k- Then, it means that we 
look at the height fluctuations at times ^(r^, 9^-, Sk) = t + Okt", with 



x(r, 0) - y{T, 9) = t2}IH^I'^, x{t, 9) + y(r, 9) = \{t + k") - {s - t^)2^^/H^/\ (12) 



Thus, one can cover much more than only the space-like regions. As before, the projec- 
tion along the characteristic line of (x, y) on 9^ = leads to f[T^ with 9 = and with 
a slightly modified r (i.e., r — )■ r(l + o(l))). Then, using slow decorrelation, one can 
extend ([6]) to the following result: fix m > 1, z/ G [0, 1), real numbers Ti < T2 < ■ ■ ■ < Tm 
and Si, . . . , Sm- Set 



(11) 



l.h.s. of dS]) 



v=0 




X{r) = t2'IH^'\ H{t, 9,s) = ^{t + 9^) + (r^ - s)2-'/H'/\ 
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Then we have 

/ m 

InnP I f]{h{X{n),t + ekt'')>H{rk,ek,Sk)} ) = P ( []{A2in) < Sj,} ] . (13) 



\k=l J \k=l 

This type of computations can be readily adapted to the other KPZ models considered 
in the sequel. 

3.2. Point to point LPP. Consider the following random measure fi on M'^"'"^: 

= ^ Wp6p 

where Z_|_ = {0, 1, . . .} and Wp are non-negative random variables. One may consider 
directed paths to be restricted to follow the lattice edges. This is just standard (point- 
to-point) last passage percolation (as considered for instance in [23] )• We will restrict 
ourselves to the case where d = 1, i.e., LPP in the 2-dimensional corner. Here we write 
Wij for weights. 

The conditions for our slow decorrelation theorem to hold amount to the existence 
of a LLN and CLT. Presently, for point to point LPP, this is only rigorously know 
for the two solvable classes of weight distributions - exponential and geometric. For 
general weight distributions, the existence of a LLN follows from superadditivity (via 
the Kingman subadditive ergodic theorem), though the exact value of the LLN is not 
known beyond the solvable cases. None the less, universality is expected at the level 
of the CLT for a very wide class of underlying weight distributions. That is to say 
that, after centering by the LLN, and under t^^^ scaling, the fluctuations of LPP should 
always be given by the -Fgue distribution. In the results we now state we will restrict 
attention to exponential weights, as geometric weights lead to analogous results. 

Define LPP with two-sided boundary conditions as the model with independent ex- 
ponential weights such that, for positive parameters vr, 77 > 0, 

exponential of rate vr, if i > 0, j = 0; 

exponential of rate rj, if i = 0, j > 0; 

exponential of rate 1, if i > 0, j > 0; 

zero, if i = 0, j = 0. 



(14) 



Recall that an exponential of rate A has mean 1/A. This class of models was introduced 
in (and for geometric weights in ^) and includes the one considered in [23j. 

A full description of the one-point fluctuation limit theorems for L'^c{tp) was con- 
jectured in [31] (see Conjecture 7.1) and a complete proof was given in [6|. These limit 
theorems show that the hypotheses for slow decorrelation are satisfied and hence The- 
orem [2]T] applies. We present an adaptation of Theorem 1.3 of [6] stated in such a way 
that Theorem 12.11 is immediately applicable. As such, we also state the outcome of 
applying Theorem 12.11 (see Figure HI for an illustration). The characteristic direction u 
as well as the exponents and limiting distributions for the fluctuations depend on the 
location of the point p as well as the values of vr and rj. Due to the radial scaling of p 
it suffices to consider p of the form p = (1, k^) for G (0, 00). 

In the following theorem, j-p-p = 1/3, c-^c is a constant depending on the direction 
p, c-p-p is a constant depending on the direction u, and D' is -Fgue- We also refer the 
reader to [6] for the definitions of the distribution functions -Fgue, Fqqe and Fq which 
arise here. 
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Theorem 3.1. Define = , = l-^ and k^h = y^^^jj^- 

(1) // > K > (which implies that tt + rj > 1) then l-^c = (1 + > luc = 1/3; 
and D is either -Fque (in the case of strict inequality) or -Fgoe T^'^ case of 
either, hut not both, equalities) or Fq (in the case where all three terms in the 
inequality are equal). Then, there is slow decorrelation in the direction u = p 
for all V G (0, 1). 

(2) If 7!' + rj > 1 and k > then i^c = ~'~ T^' '^'^^ ~ ^'^^ ^ ~ -^O-i ^'^ 
standard Gaussian distribution. Then there is slow decorrelation in the direction 
u = (1, K^) for all V e (0, 1). 

(3) If n + rj > 1 and k, < k.^ then i-nc = ^ + Inc = 1/2, and D = Aoj. Then 
there is slow decorrelation in the direction u = (1, k^) for all v G (0, 1). 

(4) If TT + ri < 1 and k > then — ^ ^ Ihc = 1/2, and D = A/o,i. Then 
there is slow decorrelation in the direction u = (1, k^) for all v G (0, 1). 

(5) If TT + rj < 1 and k < Ksh then i^c = ^ + I'HC = 1/2, and D = A/o,i. Then 
there is slow decorrelation in the direction u = (1, k^) for all v G (0, 1). 

(6) // TT + ?7 < 1 and K = then 4; = = ^^(j^, Ihi = 1/2, and D is 
distributed as the maximum of two independent Gaussian distributions. Then 
there is no slow decorrelation. 

This last passage percolation model is related to a TASEP model with two-sided initial 
conditions (which we discuss in Subsection I3.3p . As explained before the characteristics 
are those for the Burgers equation. The first three cases above correspond with a 
situation that is known of as a rarefaction fan, while the last three correspond with 
the Shockwave. The above result is illustrated in Figure ID The left case displays the 
rarefaction fan (the fanning of the characteristic lines from the origin) and the right case 
displays a shockwave (the joining together of characteristic lines coming from different 
directions). 

In addition to one-point fluctuation limits, the above two-sided boundary condition 
LPP model has a fully classified limit process description. The description was given 
in [3] for TT + rj = 1 (known as the stationary case) and in [15] for all other (non- 
equilibrium) boundary conditions. These process limits are obtained using determinan- 
tal expressions for the joint distribution of the last passage times for points along fixed 
directions. Thus, initially, one only gets process limits along fixed lines. As explained 
in Section 13.11 and in [15j slow decorrelation, however, implies that the appropriately 
rescaled fluctuations at the points which are off of this line (to order for u < 1) 
have the same joint distribution as their projection along characteristics to the line (see 
Figure 1 of |1] for an illustration of this). 

A completely analogous situation arises in the case of geometric, rather than expo- 
nential weights (this model is often called discrete PNG). Such a model is described 
in [5] and the one-point limiting fluctuation statistics are identified. The spatial pro- 
cess limit is characterized in [2T]. These results are only proved in a fixed space-time 
direction, though applying Theorem 12. II we can extend this process limit away from this 
fixed direction just as with the exponential weights. 

A slightly different model with boundary conditions was introduced in [3j and involves 
thick one-sided boundary conditions. Fix a A; G N, parameters tti, . . . ,TTk, and set tTj = 1 
for i > k. Just as above, we define independent random weights on Z^, this time 
with Wij exponential random variables of rate tTj (mean l/vTj). Section 6 of [3] explains 
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Figure 4. The different cases of fluctuation limit theorems and accom- 
panying directions of slow decorrelation: (a) 77 + vr > 1 (actually shown 
T] = IT = 2/3); (h) T] + TT < 1 (actually shown r] = tc = 1/3). As 
(the height along the dashed line) varies, the case of fluctuation theorem 
changes, as does the direction of slow decorrelation (given by the direction 
of the thin lines). 

how results they obtain for perturbed Wishart ensembles translate into a complete 
fluctuation limit theorem description for this model. Just as in the two-sided boundary 
case, those limit theorems show that the hypotheses of Theorem 12.11 are satisfied and 
therefore there is slow decorrelation. The exponent j-^c depends on the point p and 
the strength of the boundary parameters tTj and can either be 1/3, with random matrix 
type fluctuations, or 1/2 with Gaussian type (more generally i x i GUE for some 
1 < i < k) fluctuations (see |3] Theorem 1.1). The exponent •y-p-p = 1/3 and the 
limiting distribution D' is -Fgue- The direction of the slow decorrelation depends on 
the parameters and the point (we do not write out a general parametrization of this 
direction as there are many cases to consider depending on the Hi). The fluctuation 
process limit theorem has not been proved for this model, though the method of [15] 
would certainly yield such a theorem. Also, analogous results for the geometric case 
have not been proved either but should be deducible from the Schur process [H]. 

3.3. TASEP and PASEP. 

3.3.1. Totally asymmetric simple exclusion process (TASEP). TASEP is a Markov pro- 
cess in continuous time with state space {0, 1}^ (think of Is as particles and Os as holes). 
Particles jump to their right neighboring site at rate 1, provided the site is empty. The 
waiting time for a jump is exponentially distributed with mean 1 (discrete-time versions 
have geometrically distributed waiting times). See [29|[30] for a rigorous construction of 
this process. 

TASEP with different initial conditions can be readily translated into LPP with 
specific measures n and hence Theorem 12.11 may be applied. Slow decorrelation can 
thus be used to show that fluctuation limit processes can be extended from flxed space- 
time directions to general 1 + dimensional space-time directions. 

An observable of interest for TASEP is the integrated current of particles I{x,t) 
deflned as the number of particles which jumped from x to a; + 1 during the time 
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Figure 5. Rarefaction wave on the left and shockwave on the right. 

interval [0,t]. Also of interest is the height function h{x,t) 

r2/(0,t) + E-=i(l-2r/i(z)), x>l, 

= < 2/(0,t), x = 0, (15) 

[2J(0,t)-EL+i(l-2r7*«), a;<-l, 

where rit{i) = 1 (resp. rit{i) = 0) is site i is occupied (resp. empty) at time t. There is a 
simple relationship between the current and the height function given by 

I{x,t) = l{h{x,t) - x). (16) 

A well-studied initial condition is step initial condition: At time t = 0, {. . . , —2, —1, 0} 
is filled by particles and {1,2, . . .} is empty, i.e., h{x,t = 0) = There is a simple 
relation with the corner growth model studied in Subsection 13. 1[ The weight Wij is 
the waiting time (counted from the instant when the particle can jump) for the particle 
which started in position + 1 to move from position i—j to i—j + 1. Thus, the TASEP 
height function records the boundary of the region of points p for which Ly^c{p) < t. 
Therefore, as in Subsection 13. 11 one can apply Theorem 12.11 leading to slow decorrelation 
(in the sense of equation ([7j)) for the fluctuations of the TASEP height function along 
space-time lines corresponding to the characteristics of Burgers equation. 

An important generalization of the step initial condition are the two-sided Bernoulli 
initial conditions which are defined for all pairs P-, p+ € [0, 1] as the random initial con- 
ditions in which particles initially occupy sites to the left of the origin with probability 
P- (independently of other sites) and likewise to the right of the origin with probability 
p+. It was proven in ^32j that two-sided TASEP can be mappeqj to the LPP with 
two-sided boundary conditions model (|T4|) with vr = 1 — p+ and ?7 = p_. Using this 
connection and slow decorrelation, one can show that all the results stated for the LPP 
model ([T3D can be translated into their counterpart for two-sided TASEP. This is made 
in detail in [13] (which uses some arguments of this paper), where we prove a complete 
fluctuation process limit for p_ 7^ p+ which complements the recent result of [|] for 

P- = P+- 

The characteristic line leaving position x has slope 1 — 2p{x). On top of this, the 
entropy condition ensures that if p^ > p^, there will be a rarefaction fan from the origin 
which will fill the void between lines of slope 1 — 2p_ and 1 — 2p+. The Rankine-Hugoniot 
condition applies to the case where p_ < p+ and introduces Shockwaves with specified 



The mapping requires a geometric number of zero weights along the boundary which do not affect 
asymptotics. 
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(a) (b) 

Figure 6. Flat (a) and half- flat (b) TASEP correspond, via their height 
functions, with the LPP models in the two regions shown above. Char- 
acteristic lines are perpendicular to the initial direction of the height 
function and in case (b) the entropy condition implies that they fan out 
to the right of the origin. 

velocities when characteristic lines would cross. These two types of characteristics are 
illustrated side- by-side in Figure |5l 

Another variation is TASEP with slow particles or slow start-up times, which is con- 
sidered in [2]. It may likewise be connected to the LPP with thick one-sided boundary 
conditions model which we previously introduced. As a result we may similarly conclude 
slow decorrelation. 

Not all initial conditions correspond to LPP with weights restricted to Z^. For 
example, TASEP with flat (or periodic) initial conditions corresponds to the case where 
only the sites of fcZ, for k > 2 are initially occupied. For simplicity, we focus on the case 
k = 2. Then, the height function at time t = is a saw-tooth, see Figure [6]^a) (though 
asymptotically fiat, from which the name). Rotating by 7r/2, it is the growth interface 
for half-line to point LPP where the measure fi is supported on points in {i,j) G 
such that i + j > and given by delta masses with independent exponential weights of 
rate 1. Fluctuation theorems and limit process have been proved for several periodic 
initial conditions [Hllin] (in [3 was in discrete time, i.e., with geometric weights). 

Similarly TASEP with half-flat initial conditions is defined by letting particles start 
at 2Z_ = {■■■, —4, —2, 0}. The corresponding last passage percolation model has non- 
zero weights for points {i,j) such that i + j > and j > 0. The limit process for this 
model was identified in [IT]. Theorem 12 . 1 1 applies to both of these model and proves slow 
decorrelation. This implies that the fluctuation process limits extend to general 1 + 
dimensional space-time directions. The characteristics lines are shown in Figure [6]^b). 

A variant of half fiat initial conditions has particles starting at 2Z_ plus a few particles 
at positive even integers, with a different speed a. This is known as two speed TASEP 
and [13] gives a complete description and proof of the process limit for these initial 
conditions. As with all of the other examples, this can be coupled with a LPP model 
and hence Theorem 12.11 applies and prove slow decorrelation and enables us to extend 
these process limit results as well. 

3.3.2. Partially asymmetric simple exclusion process (PASEP). The PASEP is a gener- 
alization of TASEP where, particles jump to the right-neighboring site with rate p ^ 1/2 
and to the left-neighboring site with rate q = 1 — p (always provided that the destina- 
tion sites are empty). An important tool to study PASEP is the basic coupling [291130] . 
Through a graphical construction, one can realize and hence couple together every 
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PASEP (with different initial conditions) on the same probabihty space. Even though 
PASEP can not be mapped to a LPP model, it still has the same super-additivity prop- 
erties necessary to prove a version of Theorem 12.11 The property comes in the form 
of attractiveness. That PASEP is attractive means that if you start with two initial 
conditions corresponding to height functions hi{x,0) < h2{x,0) for all a; G M, then for 
any future time t, hi{x,t) < h2{x,t) for all x G M. 

We now briefly review this graphical construction. Above every integer draw a half- 
infinite time ladder. Fix p (and hence q) and for each ladder place right and left 
horizontal arrows independently at Poisson points with rates p and q respectively. This 
is the common environment in which all initial conditions may be coupled. Particles 
move upwards in time until they encounter an arrow leaving their ladder. They try to 
follow this ladder, and hence hop one step, yet this move is excluded if there is already 
another particle on the neighboring ladder. That this graphical construction leads to 
attractiveness for the PASEP is shown, for instance, in [29|l30]. 

In a series of three papers [381410] Tracy and Widom show that for step initial con- 
ditions with positive drift 7 = p — g > 0, PASEP behaves asymptotically the same as 
TASEP (when speeded- up by I/7). Just as in TASEP the current or height function is 
of central interest. I{x, t) is defined as the number of particles which jumped from x to 
X + 1 minus the ones from x + 1 to x during [0,t] and h{x,t) is defined by (1151) . This 
time, the height function does not monotonically grow, but does still have a drift. 

The slow decorrelation theorem for PASEP with general initial conditions is stated 
below. By a PASEP model we mean a measure on initial configurations, as well as a 
rate p = 1 — g G (1/2,1]. We write h{x,t) for the height function for this specified 
model, and h'{x,t) for the height function for the PASEP with step initial conditions. 
Note that the generalizations of Remark 12.21 apply in this case too. 

Theorem 3.2 (Slow decorrelation for PASEP). Consider a velocity G M and a second 
variable m G M. // there exist constants (depending on v and u and the model): I and ^ 
non-negative; 7,7' G (0, 1); v G (0,7/7'); o-^^^^ distributions D and D' such that 

h(vt,t)-t£ 

— =^ U, as t goes to infinity, 

h{vt + ut" ,t + t") - tl - tn' 



t-y 

h'{ut,t)-te' _ 



D, as t goes to infinity, (17) 



D' , as t goes to infinity. 



then we have slow decorrelation of the PASEP height function at speed v, in the direction 
given by u and with scaling exponent v, i.e., for all M > 0, 

lim W{\h{vt + ue, t + f)- h{vt, t) - > Mf^) = 0. 

Proof of Theorem \3.2 . Rather than the height function we focus on the current which 
is related via equation (ITB]) . I{vt + ut^ ,t + t^) is equal to the current I{vt,t) up to 
time t, plus the current of particles which cross the space-time line from vt at time t 
to vt + ut" at time t + t^ . We consider a coupled system starting at time t reset so as 
to appear to be in step initial conditions centered at position vt. By attractiveness of 
the basic coupling, the current across the space-time line from vt at time t to vt + ut^ 
at time t + for this "step" system will exceed that for the original system. Denote 
by I'{ut^,t^) the current associated to the coupled "step" system and observe that, it 
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is distributed as the current of an independent step initial condition PASEP. Thus, 

I{vt + ut\ t + t'') = I{vt, t) + /'(wt^ t") + Xt 

where < 0. From this point on the proof follows exactly as in the proof of Theo- 
rem O □ 

Using the fluctuation results proved in [381 - HO] . reviewed in |41j, we find that the 
above hypotheses are satisfied for PASEP with step initial conditions and also for step 
Bernoulli initial conditions [l2] with p_ = and p+ > 0. The slow decorrelation 
directions are given by the characteristics just as in the case of TASEP. These two 
sets of initial conditions are the only ones for which fluctuations theorems are presently 
known for PASEP, but limit process theorems are not yet proven. 

3.4. The polynuclear growth (PNG) model. As mentioned before, slow decorrela- 
tion for the (continuous time, Poisson point) PNG model was proved previously in [19] 
in the case of the PNG droplet and stationary PNG. Theorem 12.11 (along with the nec- 
essary preexisting fluctuation theorems) gives an alternative proof of these results as 
well as the analogous result for flat PNG. Because of the minimality of the hypotheses 
of our theorem we may further prove slow decorrelation for the model of PNG with two 
(constant) external sources considered in [5]. The way that PNG fits into the framework 
of our half-line to point LPP model is that one takes /z to be a Poisson point process of 
specified intensity on some domain. For the PNG droplet, stationary PNG and PNG 
with two external sources, we restrict the point process to and (in the second and 
third cases) augment the measure with additional one dimensional point process along 
the boundaries. For fiat PNG the support of the point process is {(x, y) : x + y > 0}. 
The limit process for the PNG droplet for fixed time was proved in [32] and for fiat 
PNG was proved in [T2] for space-like paths. It was explained in [19] that slow decorre- 
lation implies that these limit processes extend to general 1-1-0 dimensional space-time 
directions (with time scaling for z/ < 1). 

3.5. First passage percolation. As opposed to LPP one can look to the minimum 
value of T{tt). This then goes by the name of directed first passage percolation and for 
simplicity we consider this only when we restrict our measure to being supported on a 
lattice. One may also consider undirected first passage percolation. Theorem 12.11 can 
be adapted in a straightforward way for both of these models. The statement of the 
theorem remains identical up to replacing the last passage time variable with the first 
passage time. For the proof the only change is that the compensator Xt now satisfies 
Xt < rather than Xt > 0. Unfortunately no fluctuation theorems have been proved 
for first passage percolation, so all that we get is a criterion for slow decorrelation. 

3.6. Directed polymers. We now briefly consider a lattice-based directed polymer 
models in 1 -|- 1 dimension and note that just as in LPP, slow decorrelation can arise in 
these models. Unfortunately, just as in first passage percolation, there are no fluctuation 
theorems proved for such polymers. Recently, however, the order of fluctuations for a 
particular specialization of this model was proved in [35]- It should be noted that while 
we focus on just one model, the methods used can be applied to other polymer models 
and in more than 1-1-1 dimension (for example line to point polymers). 

The model we consider is the point to point directed polymer. In this model we 
consider any directed, lattice path tt from (0, 0) to a point p and assign it a Gibbs 
weight e^'^^'^^ where /3 > is known as the inverse temperature and where T(7r) is the 
sum of weights (which are independent) along the path vr (— T(7r) is the energy of the 
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path vr). We define the partition function and free energy for a polymer from a point p 
to q as: 

Zp{p,q) = Yl e^^^'^^ = bogZ^{0,p). 

It is expected that the free energy satisfies similar fluctuation theorems to those of LPP 
(which is the (3 = oo limit of 5^/3 (0,p)). 

Theorem 3.3. Consider a directed polymer model and consider a point p G and a 
direction u G . // there exist constants ( depending on p and u and the model weight 
distributions): i and i' non-negative; 7,7' G (0,1); v G (0,7/7'); ^'^^ distributions D, 
D' such that 

Xi{t) '■= — =^ D, as t goes to infinity, 

m d^{0,tp + t'^u)-U-tn' , • ^ V 

X2{t) '■= — — =^ U, as t goes to infinity, 

^p{tp,tp + t''u)-t-i' , • ^ V 

X-i[t) ■■= j-py^, ^ D , ast goes to infinity, 

then we have slow decorrelation of the point to point polymer at tp, in the direction u 
and with scaling exponent v, which is to say that for all M > 0, 

lim mdi^iO, tp + Tm) - df^{0, tp) - I > Mt'^) = 0. 

Proof. The direction u for a given p should correspond to the characteristic through 
that point. The proof of this criterion for slow decorrelation is identical to the proof for 
Theorem 12.11 and follows from the computation below (a result of super-additivity yet 
again) : 

MO, tp + fu) = ^ log ( E e^^^-^ X ^'^^^^ + E ^'^^'^O 

^ ^7r:0-S>tp -K-.tp-^tp+fu TT-.O^tp+fu, ' 

tp(fi1T 



\^og(Y ^'''^^"^^ E e^^^A+Xt 



/5 

^ 7r:0— S>tp TT-.tp-^tp+t'-'u 

= dfsiO, tp) + dp{tp, tp + fu) + Xt. 
Here > and the argument is analogous to ([3]). □ 
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